Optimal flexibility for conformational transitions in macromolecules 
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Conformational transitions in macromolecular complexes often involve the reorientation of lever- 
like structures. Using a simple theoretical model, we show that the rate of such transitions is 
drastically enhanced if the lever is bendable, e.g. at a localized "hinge" . Surprisingly, the transition 
is fastest with an intermediate flexibility of the hinge. In this intermediate regime, the transition 
rate is also least sensitive to the amount of "cargo" attached to the lever arm, which could be 
exploited by molecular motors. To explain this effect, we generalize the Kramers-Langer theory for 
multi-dimensional barrier crossing to configuration dependent mobility matrices. 



Many biological functions depend on transitions in the 
global conformation of macromolecules, and the associ- 
ated kinetic rates can be under strong evolutionary pres- 
sure. For instance, the directed motion of molecular mo- 
tors is based on power strokesll], protein binding to 
DNA can require DNA bending [2| or spontaneous partial 
unwrapping of DNA from histones [J El) an d t ne func- 
tioning of some ribozymes depends on global transitions 
in the tertiary structure [5[. These and other examples 
display two generic features: (i) A long segment within 
the molecule or complex is turned during the transition, 
e.g. an RNA stem in a ribozyme, the DNA as it unwraps 
from histones or bends upon protein binding, or the lever 
arm of a molecular motor relative to the attached head, 
(ii) The segment has a certain bending flexibility. Here, 
we use a minimal physical model to study the coupled 
dynamics of the transition and the bending fluctuations. 

Our model, illustrated in Fig. [TJ demonstrates explic- 
itly how even a small bending flexibility can drastically 
accelerate the transition. Furthermore, if the flexibil- 
ity arises through a localized "hinge", e.g. in the protein 
structure of some molecular motors [6, 7[ or an interior 
loop in an RNA stem, we find that the transition rate 
is maximal at an intermediate hinge stiffness. Thus, in 
situations where rapid transition rates are crucial, molec- 
ular evolution could tune a hinge stiffness to the optimal 
value. We find that an intermediate stiffness is optimal 
also from the perspective of robustness, since it renders 
the transition rate least sensitive to changes in the drag 
on the lever arm, incurred e.g. by different cargos trans- 
ported by a molecular motor. 

Our finding of an optimal rate is reminiscent of a phe- 
nomenon known as resonant activation 0, where a 
transition rate displays a peak as a function of the charac- 
teristic timescale of fluctuations in the potential barrier. 
However, we will see that the peak in our system has 
a different origin: a trade-off between the accelerating 
effect of the bending fluctuations and a decreasing av- 
erage mobility of the reaction coordinate. The standard 
Kramers-Langer theory [Io| for multi-dimensional transi- 
tion processes is not sufficient to capture this trade-off. A 



generalization of the theory to the case of configuration- 
dependent mobility matrices turns out to be essential to 
understand the peak at intermediate stiffness. 

Model. — We model the conformational transition as 
a thermally activated change in the attachment angle tp 
of a macromolecular lever, see Fig. [T] The lever has two 
segments connected by a hinge with stiffness e, which ren- 
ders the lever preferentially straight, but allows thermal 
fluctuations in the bending angle 9. The energy function 
V(ip, 9) of this 'Two-Segment Lever' (TSL) is 
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where ksT is the thermal energy unit. The hinge, de- 
scribed by the first term, serves not only as a sim- 
ple model for a protein or RNA hinge, but also as a 
zcroth-order approximation to a more continuously dis- 
tributed flexibility; see below. The second term is the 
potential on the attachment angle ip 1 which produces a 
metastable minimum at (<p,0) = (0,0). The thermally- 
assisted escape from this minimum passes through the 
transition state at (<p, 9) = (b/a,0) with a barrier height 
AV = b 3 k B T/6 0. 
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FIG. I: Schematic illustration of the 'Two-Segment Lever' 
(TSL) model for conformational transitions, (a) The two 
segments of lengths I and p are connected by a hinge and 
attached to the origin. The viscous drag acts on the ends 
of the segments as indicated by the beads, (b) Schematic 
illustration of the barrier crossing processes. The external 
meta-stable potential V{ip) is indicated by shading (top; dark 
corresponds to high energy) and is also sketched below. 
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In the present context, inertial forces are negligible, 
i.e. it is sufficient to consider the stochastic dynamics of 
the TSL in the overdamped limit. We localize the fric- 
tion forces to the ends of the two segments, as indicated 
by the beads in Fig.QTa). The length of the first segment 
defines our length unit and p denotes the relative length 
of the second segment. Similarly, we choose our time 
unit such that the friction coefficient of the first bead is 
unity, and denote the coefficient of the second bead by 
£. To describe the Brownian dynamics of the TSL, we 
derive the Fokker-Planck equation for the time-evolution 
of the configurational probability density p(ip,9,t). In 
general, the derivation of the correct dynamic equations 
can be a nontrivial task for stochastic systems with con- 
straints U,[3| . For instance, implementing fixed seg- 
ment lengths through the limit of stiff springs, leads to 
Fokker-Planck equations with equilibrium distributions 
that depend on the way in which the limit is taken [1^ . 
However, for our overdamped system, we can avoid this 
problem by imposing the desired equilibrium distribu- 
tion, i.e. the Boltzmann distribution p = exp(— V/k^T), 
which together with the well-defined deterministic equa- 
tions of motion uniquely determines the Fokker-Planck 
equation for the TSL. 

The deterministic equations of motion take the form 
<jk = — Mfci dV/dqi with the coordinates (31, 52) = (<P, 0) 
and a mobility matrix M. We obtain M with a standard 
Lagrange procedure: Given linear friction, M is the in- 
verse of the friction matrix, which in turn is the Hessian 
matrix of the dissipation function [l3| • This yields 
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The Fokker-Planck equation then follows from the conti- 
nuity equation dtp({qi}, t) = —dkjk({qi}, t) together with 
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as the probability flux density. Our analytical analysis 
below is based directly on Eqs. (0) and (J3]), while we 
perform all Brownian dynamics simulations with a set of 
equivalent stochastic differential equations 14 1. 

Transition rate. — To explore the phenomenology of 
the TSL, we performed simulations to determine its av- 
erage dwell time r in the metastable state, for a range of 
hinge stiffnesses e. The rate for the conformational tran- 
sition is related to the dwell time by k(e) = l/r(e). Fig. [2] 
shows k(e) (circles) for a barrier AV = 12 fcsT, a distance 
A<f = 0.4 to the transition state, and £ = p = 1 (data for 
different parameter values behaves qualitatively similar, 
as long as the process is reaction-limited, i.e. AV is suf- 
ficiently large that r is much longer than the time for the 
TSL to freely diffuse over an angle Aip). We observe a 
significant flexibility-induced enhancement of the transi- 
tion rate over a broad range of stiffnesses, compared to 
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FIG. 2: Simulation data of the barrier crossing rate normal- 
ized by ko display a prominent peak at finite stiffness (cir- 
cles, each obtained from 20000 simulation runs initialized at 
the metastable minimum). The conventional Langer theory 
fails to describe the non-monotonicity of the rate and over- 
estimates the rate at small e. The generalized Langer theory 
captures the non-monotonicity of the rate and describes the 
simulations data accurately; parameters see main text. 



the dynamics in the stiff limit (e— * 00), see inset. Note 
that the enhancement persists even at relatively large e, 
where typical thermal bending fluctuations 5<p ~ y/e are 
significantly smaller than Aip. Surprisingly, the acceler- 
ation is strongest at an intermediate stiffness (e 10). 
This observation suggests that the stiffness of molecular 
hinges could be used, by evolution or in synthetic con- 
structs, to tune and optimize reaction rates. 

When the friction coefficient £ of the outer bead is 
increased, the rate of the conformational transition de- 
creases; see Fig. [3^. This decrease is most dramatic in 
the stiff limit (dash-dotted line) . In the flexible limit (dia- 
monds) the decrease is less pronounced. Notably, the rate 
appears least sensitive to the viscous drag on the outer 
bead at intermediate e (circles). Indeed, Fig. [SJd shows 
that the e-dependence of this sensitivity (measured as the 
slope of the curves in Fig. [3^ at £= 1) has a pronounced 
minimum at e w 20. Hence, intermediate hinge stiffnesses 
in the TSL lead to maximal robustness, which is an im- 
portant design constraint for many biomolecular mecha- 
nisms in the cellular context. For instance, as molecular 
motors transport various cargos along one-dimensional 
filaments, it may be advantageous to render their speed 
insensitive to the cargo size, e.g. to avoid "traffic jams". 

In the remainder of this letter, we seek a theoretical 
understanding of the above phenomenology. First, it is 
instructive to consider simple bounds on the transition 
rate. An upper bound is obtained by completely elimi- 
nating the outer bead. The Kramers rate [l5| for the re- 
maining ID escape process, fco = (a 2 fe/27r) e~ Ay / feEiT , is 
used in Figs. and [3] to normalize the transition rates. At 
the optimal stiffness, the transition rate in Fig. [5] comes 
within 20 % of this upper bound. An obvious lower bound 
is the stiff limit: For e— >oo, the second segment increases 
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FIG. 3: The sensitivity of the rate to the friction coefficient 
£ is minimal at intermediate stiffness, (a) Simulation results 
at e = and e = 25 as well as the theoretical estimates of the 
rate at e = and in the stiff limit, (b) The derivative of In k 
with respect to ln£ evaluated at £ = 1, i.e. the slope of the 
curves in a), is minimal in an intermediate stiffness range. 




FIG. 4: The friction opposing rotation of the attachment an- 
gle tp depends on the bending angle 9, since the outer bead is 
moved by different amounts in different configurations. For an 
infinitesimal displacement dtp, the displacement of the outer 
bead is sin# dtp. The projection of the resulting friction force 
onto the direction of motion adds another factor sin 9, yielding 
a friction coefficient for ip of I + £sin 2 9. 



the rotational friction by a factor ( = 1 + (1 + p) 2 £, so 
that the ID Kramers rate becomes fcoo = ko/(, as shown 
by the dash-dotted line in Fig. [2] and Fig. [3^- However, 
to understand how the dynamics of the bending fluctua- 
tions affects the transition rate, we must consider the full 
2D dynamics of the TSL. The multi-dimensional gener- 
alization of Kramers theory is Langer's formula for the 
escape rate over a saddle in a potential landscape jlOj, 

A / det e( w ) / AV\ 
r = - x y exp \- — j . (4) 

Here, and e^ s ^ denote the Hessian matrix of the po- 
tential energy, d 2 V / dqkdqi, evaluated at the well bot- 
tom and the saddle point, respectively, whereas A is the 
unique negative eigenvalue of the product of the mobility 
matrix M and e( s ) . Eq. ((4J can be made plausible in sim- 
ple terms: Given a quasi-equilibrium in the metastable 
state, the second factor represents the probability of be- 
ing in the transition region, i.e. the region within ~ k^T 
of the saddle. The escape rate is then given by this prob- 
ability multiplied by the rate A at which the system re- 
laxes out of the transition state, analogous to Michaelis- 
Menten reaction kinetics. 

For our potential |T]), the determinants in (01 can- 
cel. The eigenvalue can be determined analytically (the 
dashed line in Fig. [2] shows the resulting fcLangcr), but for 
the present purpose it is more instructive to consider the 
expansions for large and small stiffness, fn the stiff limit, 
the natural small parameter is the stiffness ratio 7/e, 
where 7 = a 2 b is the absolute curvature or "stiffness" of 
the external potential at the transition state. The ex- 
pansion yields /CLanger/fcco = 1 + (p 2 £/C) 7/ e + 0(j 2 /e 2 ). 
As expected, the rate approaches , but the stiff limit 
is attained only when the bending fluctuations ~ y/e are 
small compared to the width of the barrier ~ ^7. fn the 
opposite limit, e <C 7, the rate is given by fcLanger/fco = 
I — (l + p -1 ) 2 e/7 + 0(e 2 /j). Since the linear term is 
negative, Langer theory predicts that the transition rate 
peaks at zero stiffness, with a peak value equal to the 
Kramers rate fco for the lever without the second segment. 



This prediction is clearly at variance with the simulation 
results, ft is interesting to note, however, that the slope 
of the linear decay is independent of £. This is consistent 
with our observation that the transition rate is insensi- 
tive to £ in the intermediate stiffness regime. Indeed, 
Fig. [5] shows that Langer theory (dashed line) describes 
the simulation data (circles) reasonably well for interme- 
diate and large hinge stiffness. 

To understand the origin of the peak at intermediate 
stiffness, it is useful to consider the flexible limit (e = 0). 
In this limit, the transition state is degenerate in 6, and 
it seems plausible to estimate the transition rate by using 
a ^-averaged mobility for the reaction coordinate tp, 

^=0)^0/"^=^=. ( 5 ) 

This estimate agrees well with the simulation data, 
see the dashed line in Fig. [3^, indicating that the 
configuration-dependent mobility |2]) plays an important 
role for the transition rate. In contrast, the conventional 
Langer theory assumes the mobility matrix to be con- 
stant in the relevant region near the transition state. 
Fig. 0] illustrates why the mobility M u of the coordinate 
tp is affected by the bending angle 6 and gives a graphical 
construction for Mix- 
Generalized Langer theory. — To account for the mo- 
bility effect identified above, we must generalize the 
Langer theory to configuration-dependent mobility ma- 
trices. The special case where the mobility varies only 
along the reaction coordinate has already been studied in 
1 1 61 ] , however the main effect in our case is due to the vari- 
ation in the transverse direction. In the following, we out- 
line the derivation of the central result, while all details 
will be presented elsewhere. Near the saddle, the mobil- 
ity matrix takes the form Mij({qi}) = M ( y + ^A^qiq k , 
where qi are deviations from the saddle and denotes 
the tensor of second derivatives of the mobility matrix 
(we assume that the first derivatives of M vanish at the 
saddle, which is the case for the TSL). The escape rate 
is given by the probability flux out of the metastable 
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well, divided by the population inside the well. To cal- 
culate the flux, we construct a steady state solution to 
the Fokker-Planck equation in the vicinity of the saddle, 
as described in [l5| for the conventional Langer theory. 
We use the Ansatz p({qi}) = \peq{{<li}) erfc(u), where 
PeqiUi}) = Z- 1 e- 1/({ ' ?i})/fcBT and erfc(u) is the comple- 
mentary error function with argument u = U^qk- Insert- 
ing the Ansatz into the Fokker-Planck equation yields an 
equation for the vector U, 



Uii-Mijeft + B lk ) - UtMijUj U k = 



(6) 



where B ik = k^T^2, n A^. B^qk is the noise induced 
drift, which is absent in the conventional Langer theory. 
Ignoring higher order terms, this equation determines U 
to be the left eigenvector of — M( S W S ' + B to the unique 
positive eigenvalue A, and requires U to be normalized 
rWi 



such that UiMf/Uj = A 



The directions of the left and 
right eigenvectors of — M^ s ^e^ + B have a physical in- 
terpretation: U is perpendicular to the stochastic sepa- 
ratrix, while the corresponding right eigenvector points 
in the direction of the diffusive flux at the saddle [17j • 

From p({qi}), the flux density is determined by |3]) and 
the total flux is obtained by integrating the flux density 
over a plane containing the saddle; a convenient choice is 
the plane u = 0. Evaluation of the integral is particularly 
simple in a coordinate system, where the first coordinate 
is parallel to U, and the remaining coordinates are chosen 
such that e( s ) is diagonal in this subspace, = fXiSij 
for i, j > 1. In this coordinate system, the generalized 
Langer rate takes the simple form 



k = 



\ i j L_ V ^JJ- 



2tt 



where c = Ui&y Uj + 1 = 
the inverse matrix of e' s ) 



det e( w ) 
I det e( s ) I 



Bue^ 1 /M[f and e 1 denotes 
Eq. ([7]) contains three correc- 
tions to (jlj, all of which vanish when M({qi}) is constant: 
The most important one is given by X)z>i ^ii/W> which 
changes the mobility M\\ in the direction of U to an ef- 
fective mobility that is averaged over the separatrix with 
respect to the Boltzmann distribution. In addition, there 
are two corrections incurred by the noise-induced drift: 
the factor y/1 — c and a change due to the fact that A is 
now the eigenvalue to M( S W S ) — B instead of M^e^. 

The solid line in Fig. [2] shows the application of the 
generalized Langer formula to the TSL. We observe that 
it captures the peak in the transition rate and thus the es- 
sential phenomenology of the TSL. Obviously, the evalu- 
ation of the Gaussian integral that leads to Eq. |(7J) is only 
meaningful, if the harmonic approximation of the mobil- 
ity matrix is reasonable within the relevant saddle point 
region. This integral diverges as the saddle point degen- 
erates, which explains the behavior for e — » 0. At high 
£, the very anisotropic friction can also render Langer 
theory invalid HI, 13 |. 



(7) 



Conclusion. — We have introduced the "Two-Segment 
Lever" as a simple model for a class of conformational 
transitions in biomolecules. The model clearly demon- 
strates how flexibility can enhance the rate of a confor- 
mational transition. This remains true, if the hinge in the 
TSL is replaced by a more continuous bendability. In- 
terestingly, a discrete hinge has a stiffness regime, where 
the rate is large and robust against cargo variation, which 
raises the question, whether these effects are exploited by 
evolution, for example in the design of molecular motors. 
To understand these effects theoretically, we derived a 
generalized Langer theory that takes into account con- 
figuration dependent mobility matrices. We hope that 
this theory will find applications also in other fields. 
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